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Abstract. Geometric generalized Mittag-Leffler distributions having the Laplace transform 



l 



,0<a<2,/3>0is introduced and its properties are discussed. Autoregressive 



CN . l+/31og(l+t«)> 

processes with Mittag-Leffler and geometric generalized Mittag-Leffler marginal distributions are 
developed. Haubold and Mathai (2000) derived a closed form representation of the fractional ki- 
netic equation and thermonuclear function in terms of Mittag-Leffler function. Saxena et al (2002, 
2004a, b) extended the result and derived the solutions of a number of fractional kinetic equations 
in terms of generalized Mittag-Leffler functions. These results are useful in explaining various fun- 
damental laws of physics. Here we develop first-order autoregressive time series models and the 
properties are explored. The results have applications in various areas like astrophysics, space sci- 
ences, meteorology, financial modeling and reliability modeling. 

Keywords: Autoregressive process, o-Laplace distribution, geometric infinite divisibility, geometric generalized 
Mittag-Leffler distribution, generalized Mittag-Leffler distribution, self decomposability, time series modeling, finan- 
cial modeling. 

a 

1 Introduction 

' Recently Mittag-Leffler functions and distributions have received the attention of mathematicians, 

statisticians and scientists in physical and chemical sciences. Pillai (f990) introduced the Mittag- 
| Lefflcr distribution in terms of Mittag-Leffler functions. Jayakumar and Pillai (1993) developed a 

first order autoregressive process with Mittag-Leffler marginal distribution. Fujita (1993) discussed 
a generalization of the results of Pillai. Jose and Seetha Lekshmi (1997) developed geometric ex- 
ponential distribution and Seetha Lekshmi and Jose (2002, 2003) extended the results to obtain 
geometric Mittag-Leffler distributions. Jayakumar and Ajitha (2003) obtained various results on ge- 
ometric Mittag-Leffler distributions. Seetha Lekshmi and Jose (2004) extend the concept to obtain 
geometric a-Laplace processes. 

A discrete version of the Mittag-Leffler distribution was introduced by Pillai and Jayakumar 
(1995). Lin (1998a, b, 2001) obtained various characterizations of the Mittag-Leffler distributions 
and refers to them as positive Linnik laws. In physics, Haubold and Mathai (2000) derived a closed 
form representation of the fractional kinetic equation and thermonuclear function in terms of Mittag- 
Leffler function. Saxena et al. (2004a, b) extended the result and derived the solutions of a number 
of fractional kinetic equations in terms of generalized Mittag-Leffler functions and obtained the 
solution of a unified form of generalized fractional kinetic equations, which provides the unification 
and extension of the earlier results. Such behaviors occur frequently in chemistry, thermodynamical 
and statistical analysis. In all such situations the solutions can be expressed in terms of generalized 
Mittag-Leffler functions. Weron and Kotulski (1986) use Mittag-Leffler distribution in explaining 
Cole-Cole relaxation. 

was first introduced by Mittag-Leffler in 1903 (Erdelyi, 



The function E a (z) — 

fc=0 
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_T(l + ak) 

1955). Many properties of the function follow from Mittag-Leffler integral representation 
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1 f t a ^ 1 e t 
EM) = 7— / dt. 

where the path of integration C is a loop which starts and ends at — oo and encircles the circular 
disc \t\ < z~. Pillai(1990) proved that F a (x) = 1 - E a (-x a ),0 < a < 1 are distribution functions, 
having the Laplace transform ip(t) = (1 +t a )~ 1 ,t > which is completely monotone for < a < 1. 
He called F a (x), for < a < 1, a Mittag-LefHer distribution. The Mittag-Leffler distribution is 
a generalization of the exponential distribution, since for a = 1, we get exponential distribution. 
Pillai (1990) has shown that F a (x) is geometrically infinitely divisible (g.i.d.) and is in the domain 
of attraction of stable laws. 

Pillai (1985) developed a-Laplace distribution with characteristic function given by (1 + |t|") _1 ; 
< a < 2. This distribution is also known as Linnik distribution. Using a simple characterization 
of this distribution given by Linnik (1962), Anderson and Arnold (1993) constructed a discrete- 
time process having a stationary Linnik distribution. These Linnik models appear to be viable 
alternatives to stable processes as models for temporal changes in stock prices. Pakes (1998) gave 
a mixture representation for symmetric generalized Linnik laws as A a _ v = S a , where A a , v 

denotes a random variable having the characteristic function (1 + \t\ a )~ v ,v > 0. This defines the 
generalized symmetric Linnik law SLi(a, v). Kotz and Ostrowski(1996) proved the following mixture 
representation of Linnik law in terms of another: If < j3 < a < 2, then there is a random variable 
V a ,p{v) > whose density is 

a . 7T/3 x' 3 - 1 
g(x,a,P) = -sin( — )- 



7T a J l + x 2 f3 + 2xf 3 cos{^-) 



such that Y = XV a ,p{v) where X = AL(a); Y = AL(f3). Jacques et al (1999) proved that the 
generalized Linnik laws belong to the Paretian family. They have shown that 

V 7TCt 

lim x a P[X > x] = — T(a) sin( — ), where X has the Linnik distribution. They also discussed the 

ro— >oo 7T 2 

estimation of parameters of Linnik distribution. Kozubowski (2000) discussed the fractional moment 
estimation of Linnik parameters. Jayakumar et al. (1995) generalized the Laplace processes of 
Lawrance (1978) and Dewald and Lewis (1985). 

Gaver and Lewis (1980) derived the exponential solution of first order autoregressive equation 
X n = pX n -i + f] n ,n = 0, ±1,±2, where {i] n } is a sequence of independently and identically 
distributed random variables when < p < 1. Lin (1994) proved the characterizations of the 
Laplace and related distributions via geometric compounding. It is well known that the concept 
of geometric compounding is related to the rarefaction of renewal processes (Renyi, 1956) and to 
damage models (Rao and Rubin, 1964; see also Galambos and Kotz, 1978, p. 95). 

In Section 2 of this paper, we describe some of the properties of generalized Mittag-LefHer distri- 
butions. In Section 3, we develop an autoregressive process of order one (AR (1)) with generalized 
Mittag-LefHer distribution. In Section 4, geometric generalized Mittag-LefHer distribution is in- 
troduced and their properties are discussed. Section 5 deals with autoregressive processes with 
the above marginal distribution. The extension of these results to k th order case is considered in 
Section 6. In Section 7 we describe some applications in various fields. 



2 Generalized Mittag-Leffler Distribution 

In this Section we introduce a new class of distributions called generalized Mittag-Leffler distribution 
denoted by GMLD {a, ft). 

A random variable with support over (0, oo) is said to follow the generalized Mittag-Leffler distri- 
bution with parameters a and /3 if its Laplace transform is given by 
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iP(t) = E[e- tx ] = (1 + t a )- fi - < a < 1,8 > 0. (1) 
The cumulative distribution function (c.d.f.) corresponding to (1) is given by 

^ (-l) k T(B + k)x a( - f3+k ^ 

It easily follows that when 8=1, we get Pillai's Mittag-LefHer distribution (see Pillai, 1990). When 
a=l, we get the gamma distribution. When a=l, B—l we get the exponential distribution. This 
family may be regarded as the positive counterpart of Pakes generalized Linnik distribution charac- 
terized by the characteristic function 

(l + |t| Q )^;0<a<2,/3>0. (3) 

(see Pakes, 1998). Now we shall discuss some properties of generalized Mittag-LefHer distributions. 

Theorem 2.1 

Let U a follows the positive stable distribution with Laplace transform ij)(t) = exp(—t a ), 
t>0;0<a<l. Let Vp be a random variable, independent of U a , and follows a gamma distribution 

/ N/3 i 

with parameter 8 and Laplace transform 4>(t) = ( j^j ;8>0. Then X ai p — U a Vpo< follows 
generalized Mittag-LefHer distribution GMLD(a, ft). 

Proof: 

The Laplace transform of X is 



<Mt) = E(e~ tvUa ) 

<o 

iPu(tvi)dF(v) 

= / e- tav dF(v) 
Jo 

1 



f 

Jo 



l + t° 



Theorem 2.2 

The p.d.f. of X a .j) is a mixture of gamma densities. 
Proof: 

Using Theorem 2.1, we have the c.d.f. of X a> p as 



F a ,p{x) = / S a , v (x)dF p (v) (4) 
Jo 

where S ayV (x) is the c.d.f. of a distribution with Laplace transform exp(—vt a ) and Fp(v) is the c.d.f. 
of gamma distribution with p.d.f. f(y) = Yjjy l3 ~ 1 e~ v . 
We can rewrite (4) as 

/>oo 

F a , p (x) = / Gpix/yFdSadv) 
Jo 

Hence the p.d.f. of X is 



I 



00 a x^- 1 
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This shows that f a< p is a mixture of generalized gamma densities. For (3 — 1, f a ^ reduces to a 
mixture of Weibull densities. 

Remark 2.1 

Lin (1998) has shown that F a ^(x) is slowly varying at infinity, for a G (0, 1] and (3 > 0. 
Remark 2.2 

Pillai (1990) had obtained the fractional moments for (3 = 1, as 

E K) = ffi^El±lM;0 <,<*<!. 

In a similar manner Lin (1998) obtained the fractional moments for X a _p for < a < 1 and (3 > 

as 

cv^r x r(l - r/a)r(/3 + r/a) 

E i x ^) = r(i-r)r/3 ;-^<^< a 

= oo if r < —a/3 or r > a 



Definition 2.1: A probability distribution on R = (0, oo) is said to be in class L if its Laplace 
transform tp(t) satisfies 

i/>(t) = il>(kt)il> k (t),teR,ke (0,1), (5) 

with ipk a Laplace transform. 
Theorem 2.3 

The GMLD belongs to the class L. 
Proof: 

The proof follows because of the relation 



(1 + t a )- !) = (1 + k a t a )- fi \ k a + (1 - k a ) 




Theorem 2.4 

The GMLD (a, (5) is geometrically infinitely divisible for0<a<l,0</3<l. 
Proof: 

From Pillai and Sandhya (1990), a distribution is g.i.d. if and only if its Laplace transform is of the 
form ip(t) = 1+ 0( t ) where ip(t) has complete monotone derivative (c.m.d.) and (f)(0) = 0. Now for 
GMLD (a, f3), the Laplace transform is ip(t) = where (j)(t) = (1 + f a Y — 1- This has c.m.d. if 

and only if < a < 1 and < /? < 1. 

Remark 2.3 

Also being a mixture of gamma random variables, it is g.i.d. By Pillai and Sandhya (1990), GMLD 
is a distribution with c.m.d. Hence it is infinitely divisible. 



3 First Order Autoregressive Processes with 
GMLD(a, 13) Marginals 

Now we shall construct a first order time series model with GMLD marginals. The generalized 
Mittag-Leffler first order auto-regressive process GMLAR(l) is constituted by {X n ;n > 1} where 
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X n satisfies the equation 



X n = al„_i + r] n ; < a < 1, (7) 

where {r] n } is a sequence of independently and identically distributed random variables such that X n 
is stationary Markovian with generalized Mittag-LcfBer marginal distribution. In terms of Laplace 
transform, eq. (2) can be given as 

V»JC B (t)=K(*)^»-i(°*) ( 8 ) 

Assuming stationarity, we have, 



ipx(at) 
{l + a a t a f 
{l + t a f 

l + (at) c 
l + t a 

a a + (1 - a a ) 



1 



l + t a 



(9) 



We can regard the innovations {r] n } as the /3-fold convolutions of random variables f/„'s such that 

f with probability a Q 
" \ M„ with probability 1 — a" 

where M„'s are independently and identically distributed Mittag-LcfHcr random variables. Mittag- 
Leffler random variables can be generated easily using the following result. Let E be distributed 
as exponential with unit mean and let Q be distributed as positive stable with Laplace transform 
e~*° then X — QE~ will be distributed as Mittag-Leffler with Laplace transform (1 + t a )~ x (see 
Kozubowski and Rachev, 1999). 

Jayakumar et al. (1995) developed an algorithm to generate Linnik random variables. In a sim- 
ilar manner, the GMLAR(l) process can be generated using computers. If X = GMLAD(a, ft), 
then the process is strictly stationary. It is sufficient to verify that X n = GMLAD(a, (3) for every n. 
An inductive argument can be presented as follows. Suppose X n -i = GMLAD{a, j3). Then from 
(8) we have, 



~l + (at) a ~ 





1 




1 


l + t a 




[l + (at) a _ 




l + t a _ 



Hence the process is strictly stationary and Markovian provided X a is distributed as GMLD. 
Remark 3.1 

If Xq is distributed arbitrarily, then also the process is asymptotically Markovian with generalized 
Mittag-Leffler marginal distribution. 

Proof: 



X n — aX n _i + rj n 

n-l 

= a n X Q + ^2a k i] n - k . 
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Writing in terms of Laplace transform, 

ra-l 

V>x„W = i/>x (a n t)]l%(a k t) 



n-l r . . / fe+i .ot^P 



1 + (a K+i t) 



1 + (a k t) c 



(1 + t a y 13 as n -> oo. 



fe=0 

Hence it follows that even if X is arbitrarily distributed, the process is asymptotically stationary 
Markovian with generalized Mittag-Leffler marginals. We therefore have the following theorem. 

Theorem 3.1 

The first order autoregressive process X n = aX n _\ + rj n ,a € (0, 1) is strictly stationary Markovian 
with generalized Mittag-Leffler marginal distribution as in (1) if and only if the {rj n } are distributed 
independently and identically as the /3-fold convolution of the random variable {U n } where 



with probability a a 
M n with probability 1 — a° 



where {M n } are independently and identically distributed Mittag-Leffler random variables provided 
X = GMLD(a, (3) and independent of r\ n . 

Remark 3.2 

The model is defined for all values of 'a' such that a £ (0,1). The autocorrelation is given by 
p(r) = Cor{X n ,X n _ r ) = aW; r = 0, ±1, ±2, . . . 

3.1 Distribution of sums and bivariate distribution of (X n ,X n+ i) 
We have, 

X n+j = a 3 X n + a J ~V l+ i + a^ 2 T] n+2 + ... + r}n+j\j = 0, 1, 2, . . . 



Hence 



T r — X n + X n+ i + . . . + X„ +r _i 

r-1 

= ^[a J X„ + a J_1 77„ + i + a^ 2 r] n+2 + ... + T] n+J ] 

3=0 

r-1 



l-a r 
1 - a 



1 - a r " J 
1 -a 



Therefore the distribution of the sums T r is uniquely determined by the Laplace transform 

r-1 



i+ 



mi 



3 = 1 

r-1 

n 



1 - a r -'J 



a a + (l-a) a 



1 



l-a r -J 



The distribution of T r can be obtained by inverting the above expression. Next, the joint distribution 
of contiguous observations (X n ,X n+ i) can be given in terms of bivariate Laplace transform as, 



ipx n ,x n+1 {h,t 2 ) = E[exp(itiX n + it 2 X n+1 )} 
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= E[exp(itiX n + it 2 (aX n + /?„))] 
= E[exp(i(ti + at 2 )X n + it 2 r]n+i)] 
= ip n „(t 2 )ipx„{ti + at 2 ) 



~l + (at 2 ) a ~ 




1 





[ i + q _ 




1 + (h +at 2 ) a _ 





Since this expression is not symmetric in t\ and t 2 , it follows that the GMLAR(l) process is not 
time reversible. 



4 Geometric Generalized Mittag-Leffler distribution 

Geometric generalized Mittag-Leffler distribution is introduced and some of its properties are studied. 

Definition 4.1: A random variable Xoni?= (0, oo) is said to follow geometric generalized Mittag- 
Leffler distribution and write X = GGMLD(a, /3) if it has the Laplace transform 

Mt) = —4- -,0 < a < 2,(3 > 0. (10) 

^ w 1 + /31og(l + i a ) ~ 

Remark 4.1 

Geometric Generalized Mittag-Leffler distribution is geometrically infinitely divisible. 
Theorem 4.1 

Let X\ , X 2 , . . . are independently and identically distributed Mittag-Leffler random variables and 
N(p) be geometric with mean i, P[N(p) = k]=p(l — p) k ~ x , k = 1,2, ...,0 < p < 1. Define 

Y = X l + X 2 + . . . + X N{p) , then Y = GGMLD(a,/3). 
Proof: 

The Laplace transform of Y is 

oo 

fe=l 

1 

l + |log(l+^)' 

Hence Y = GGMLD(a, ±). 
Theorem 4.2 

Geometric generalized Mittag-Leffler distribution is the limit of geometric sum of GML(a, £) random 
variables. 

Proof: 

(1 + t a )~P = {1 + (1 + t a )~ — 1} _ " is the Laplace transform of a probability distribution since 
generalized Mittag-Leffler distribution is infinitely divisible. Hence by lemma 3.2 of Pillai (1990), 

^„(t) = {l + n[(l + t a )" -!]}"" 

, is the Laplace transform of a geometric sum of independently and identically distribute generalized 
Mittag-Leffler random variables. Taking limit as n — > oo 
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^(t) = lim ^„(t) 

= {1+ lim - I]}" 1 

n— >oo 

= [i + /3iog(i + i a )] _1 - 

Theorem 4.3 

If W and V are independent random variables such that W has geometric gamma distribution with 
Laplace transform 1+/3 \Qg^ 1+t ^ ana - V has a positive stable distribution having Laplace transform 

e"*", then W«V = U where U = GGMLD(a,/3). 
Proof: 

The Laplace transform of U is 

Mt) = E(e~ twiv ) 

i/j v (tW«)dF(w) 



POO 

L 



f , -II 

JO 



e~ wt dF(w) 



lo l + /31og(l + i«) 

5 Geometric Generalized Mittag-Leffler Processes 

In this section, we develop a first order new autoregressive process with geometric generalized Mittag- 
Leffler marginals. 

Consider an autoregressive structure given by, 

x = ( Vn with probability p 

n \ X n -i+r) n with probability (1 — p) 

where < p < 1. Now we shall construct an AR (1) process with stationary marginal distribution 
as geometric generalized Mittag-Leffler distribution GGMLD(a, (3). 

Theorem 5.1 

Consider a stationary autoregressive process {X n } with structure given by (11). A necessary and 
sufficient condition that {X n } is stationary Markovian with geometric generalized Mittag-Leffler 
marginal distribution is that {??„} is distributed as geometric Mittag-Leffler provided X is dis- 
tributed as geometric generalized Mittag-Leffler. 

Proof: 

Let us denote the Laplace transform of X n by ipx n (t) and that of r) n by ip Vn (t), equation (11) in 
terms of Laplace transform becomes 

Vx n (t) = (t) + (l - (t) 

. On assuming stationarity, it reduces to the form 

feW = P^(t) + (l-p)V*(tMj(*)- 

Writing 

^ X{t) = l + fflog(l + t") md S ° lving We g6t ' 

= 



l + /3p log(l + t°) 
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Hence it follows that r] n — GGMLD(a,p(3). 

The converse can be proved by the method of mathematical induction as follows. Now assume that 
X n _i = GGMLD(a, (3). Then 

^x n .M = ^(t)[p+(l-p)^jc n _ a (*)] 
1 



l+p/31og(l + i a ) 

= [l+^ioga+or 1 - 



p+(i-p) 



i 



l + /31og(l+t«) 



Remark 5.1 

Note that X n and ?7„ belongs to the same family of distributions. 

5.1 The joint distribution of X n and X n _i 

Consider the autoregressive structure given in (11). It can be rewritten as 

X n = I n X n -i + rj n , 
where 

P[I n = 0] = l-P[I n = l]=p,0<p<l. 
Then the joint Laplace transform of (X n ,X n -i) is given by, 

= £ ( e (iii+.W„)X n -i)^( t2 ) 

I p 1 — p 



l + /?plog(l+i£) 



l + /Hog(l + t«) l + /31og(l + ti+f 2 Q ) 



This shows that the process is not time reversible. 



6 Generalization to a k th order geometric generalized Mittag- 
Leffler Autoregressive processes 

Lawrance and Lewis (1982) constructed higher order analogs of the autoregressive eq. (11) with 
structure as given below, 



X n 



rj n with probability p 

X n -i + Vn with probability p\ 

X n -2 + Vn with probability p 2 

X n -k + Vn with probability pk 



(12) 



where pi+p2 + - • -+Pk = 1— P, < Pi,p < 1,1 = 1,2, ... ,k and r}„ is independent of {X n , X n -i, . . .}. 
In terms of Laplace transform, eq. (7) can be given as 



1pX n (t) = pip Vn (t) + Pl^X n _! {t)lp Vn (t) + P2^X„_ 2 {t)lp Vn (t) 

Assuming stationarity, we get, 

K(*) 



■ Pk^x n _ k (t)ip Vn {t) 



i>x{t) 



P +(l-p)^x(t) 

. This establishes that the results developed in Section 5 are valid in this case also. This gives rise 
to the k th order geometric generalized Mittag-Leffler autoregressive processes. 
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7 Applications 



In thermodynamical or statistical applications, one is interested in mean values of a quantity Z(t). 
Tsallis (1988) generalized the entropic functional of Boltzmann-Gibbs statistical mechanics that 
leads to q-exponential distributions. He used the mathematical simplicity of kinetic-type equations 
to emphasize the natural outcome of this distribution that corresponds exactly to the solution of 
the kinetic equation of non-linear type; the solution has power-law behavior. Saxena et al. (2004) 
showed that the fractional generalization of the linear kinetic-type equation also leads to power-law 
behavior. In both cases, solutions can be expressed in terms of generalized Mittag-Leffler functions. 
Mittag-Leffler distributions can also be used as waiting-time distributions as well as first-passage time 
distributions for certain renewal processes. Pillai (1990) developed renewal processes with geometric 
exponential as waiting time distribution. In a similar manner renewal processes with generalized 
Mittag-Leffler and geometric generalized Mittag-Leffler waiting times can be constructed. 
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